Density of kinks after a sudden quench in the quantum Ising spin chain 
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We investigate the time evolution ol the density of kinks in the spin-1/2 quantum Ising spin chain 
after a sudden quench in the transverse field strength, and find that it relaxes to a value which 
depends on the initial and the final values of the transverse field, with an oscillating power-law decay. 
We provide analytical estimates of the long-time behavior and of the asymptotic value reached after 
complete relaxation, and discuss the role of quantum criticality in the quench dynamics. We show 
that, for a dynamics at the critical point, the residual density of kinks after the quench can be 
described by equilibrium statistical mechanics at a finite temperature dictated by the energy of the 
state after the quench. On the other hand, outside of criticality it does not exhibit thermalization. 

PACS numbers: 75.40.Gb, 75.40.Cx, 75.10.Pq, 73.43.Nq 



I. INTRODUCTION 

The experimental developments in the physics of ultra- 
cold atomic gases that have been put forward in the last 
decade opened up the possibility to test some fundamen- 
tal aspects of strongly interacting quantum many-body 
systems 1 -! 2 -. In particular, by loading cold atoms in optical 
lattices, the implementation of Hubbard-like Hamiltoni- 
ans 3 - or artificial spin chain models 4 ^ has become possi- 
ble. The key features of these setups are a detailed mi- 
croscopic knowledge of the Hamiltonian, together with an 
extremely high accuracy in controlling the system param- 
eters with tunable external fields. Moreover, due to the 
surprisingly long coherence times and to small thermal 
fluctuations, the dynamics can be monitored with very 
low dissipation and noise effects 2 -. Recently emphasis has 
been put in the investigation of out-of-equilibrium prop- 
erties of strongly correlated systems. We quote for exam- 
ple the spectacular observation of a superfluid-to-Mott 
insulator quantum phase transition in optical lattices^; 
the formation of topological defects during a quench of 
trapped atomic gases through a critical point 7 -^; the ab- 
sence of thermalization in the non-equilibrium evolution 
of an integrable Bose gas^. 

Such enormous experimental potentialities raised 
an increasing theoretical interest in the study of 
the effects induced by a sudden perturbation in 
such strongly correlated systems, which goes be- 
yond a mere debate of statistical physics princi- 
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tegrability is believed to play a crucial role in the re- 
laxation to the steady state, even if a comprehensive sce- 
nario is still lacking. The nonequilibrium dynamics of 
a nonintegrable system is expected to thermalize at the 
level of individual eigenstates^Z, following the standard 
statistical mechanics prescription a 15 i 16 i 17 i 28 . On the con- 
trary, for integrable models with nontrivial integrals of 
motion the asymptotic equilibrium states usually carry 
some memory of the initial conditions and are not the 



thermal onesii^^^^^^^ 2 ^ 2 ^. In this case it 
is however possible to derive a statistical prediction for 
the steady state in terms of a generalized Gibbs ensem- 
ble 1 ^. Several works tested the relevance of such distri- 
bution for different classes of observable o 12 i 13 ' 14 i 22 ' 23 , and 
found that it indeed works, provided the observables are 
sufficiently uncorrelated with respect to the constants of 
motio n 20 ' 21 . 

In the present study we consider a quantum quench of 
the one-dimensional transverse field Ising chain; we inves- 
tigate the time evolution of the density of kinks that are 
generated after a quench in the field strength, focusing 
on the thermalization properties of this observable and 
on the role of the quantum phase transition in the quench 
dynamics. Due to the complete integrability of the Ising 
model, observables are in general not expected to follow 
the prescriptions of statistical mechanics. As a matter of 
fact, the sensitivity to integrability can be traced back 
to the property of locality of operators with respect to 
the fermion quasiparticles that diagonalize the model in 
the continuum limit. The Ising chain possesses two sec- 
tors of operators: one that is local with respect to those 
particles, and another one that is non-local 30 . Observ- 
ables belonging to the non-local sector, as the two-point 
correlation function of the order parameter, behave ther- 
mally 2 ^; on the other hand, local observables in general 
do not. Here indeed we find that the converged value 
of the density of kinks after the quench is not equal to a 
thermal expectation value, unless the system is quenched 
towards its critical point. 

The paper is organized as follows. In Sec. |TT] we de- 
scribe our model and review some standard techniques 
that are used to approach the quenched dynamics of the 
Ising chain. In Scc. lIIII wc define the density of kinks, that 
is the quantity that we are going to analyze, and give the 
prescription to evaluate it within the free fermion for- 
malism. We first focus on its behavior for the system 
at equilibrium, both at zero and at finite temperatures 
fSec. IIVI) . and then concentrate on the quenched case at 
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zero temperature (Sec.fV]), where we analyze the asymp- 
totic value as well as the finite-time transient. In Sec. IVII 
we explicitly show that the density of kinks does not ex- 
hibit a thermal behavior outside criticality, while it ther- 
malizes only at the critical point. Finally in Sec. I VIII we 
draw our conclusions. In the appendix we derive explicit 
analytic expressions for the density of kinks of the ground 
state (App. fA} and after a quench (App. [B|) . 



Statics 



The Ising Hamiltonian (fl]) can be exactly diagonalized 
by means of a Jordan- Wigner transformation (JWT), fol- 
lowed by a Bogoliubov rotatio n 32 i 33 . First one represents 
the spins in terms of JW fermions, which are defined by: 



<j 1 exp 



3 = 1 



(3) 



II. MODEL 

We consider the simplest nontrivial example of ex- 
actly solvable one-dimensional quantum many-body sys- 
tems, exhibiting a quantum phase transition, that is the 
spin-1/2 quantum Ising chain with ferromagnetic inter- 
actions 31 . This is characterized by the Hamiltonian 



w(r) = -j 



3 



[< J 3 <J 



3 + 1 + T(J l 



(1) 



where er" (a — x,y,z) are the Pauli matrices for the 
jth spin, while the parameters J and T respectively 
denote the nearest-neighbor antiferromagnetic exchange 
coupling and the transverse field strength (hereafter we 
set J = 1 as the system's energy scale and take T > 
without loss of generality; we also use units of h — ks = 
1). In this paper we will assume periodic boundary condi- 
tions for the Ising chain, and then take the limit L — > oo, 
in order to approach the thermodynamic limit. We will 
therefore suppose that the sum in Eq. (P) goes from 1 to 
L, with the rule cr£ +1 = erf and with L being the num- 
ber of spins in the chain. We assume an even number 
for L. At zero temperature, the system is a quantum 
paramagnet for T > 1 and a ferromagnet with respect to 
the coupling direction for T < 1; these two phases are 
separated by a quantum critical point at T c = 1. 

We choose to drive the system ([!} out of equilibrium 
by performing a sudden quench of the Hamiltonian pa- 
rameter T: we assume that the transverse magnetic field 
strength is suddenly quenched from Tq to T at time 
t = 0. The system is also supposed to be initialized 
in the ground state |\&(ro)) of H(T ). Due to the sudden 
variation of the field intensity, at t > the state |\l/(ro)) 
will evolve according to the new Hamiltonian Ti.(T): 



\m) = U(t) \m) = e- m(T)t l*(To)> 



(2) 



with aj — (<7j ± ierJ)/2 raising and lowering spin opera- 
tors. This leads to a Hamiltonian that is quadratic in the 
c-fcrmions. Since the ground state has an even number of 
c-fermions, and their parity is conserved during the evo- 
lution, we will assume that the actual system state \tp{t)} 
always lies in the even c-fcrmionic Hilbert space sector. 
This requires imposing antiperiodic boundary conditions 
for the fermions. 

The Hamiltonian can then be easily handled by switch- 
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mg to momentum representation: cj, — —j- c <-j , 
where the possible values of k are fixed by the an- 
tiperiodic boundary conditions, and are given by k = 
_j_ 7r(2n+l) ^ w ^ n ^ = 0, — 1. Indeed, Ji decouples 

into a sum of independent terms, each of them acting 
on the subspace (—k,+k). This is finally diagonalized 
through a Bogoliubov transformation, by introducing the 
new fermionic operators 



7i fe 



where the coefficients 



-v k c k 



-r t 

«fc C -k 



(4) 



ib k 



(•5) 



with a£ 



— 2(cos/c + r) and bk 
dispersion is characterized by 



2 sink, while their 



( a D 2 + ( b l) 2 = 2v / l + r 2 + 2rcosfc. 



(6) 



The Ising model is then recasted into a free fermionic 
system, where the Hamiltonian |T]) can be written as: 



«( r )= E 4U^+il k j-k-i) 



(7) 



0<fc<ir 



Due to the complete integrability of the Ising model, it is 
generally possible to reduce the computational cost for 
evaluating directly measurable quantities (such as dy- 
namical correlation functions or few-body observables) 
to a linear increase with the system size L. In some cases 
one can also derive an explicit analytic expression, as we 
shall explain below. 



The ground state |^(r)) of W(T) is therefore the vacuum 
state for the fermionic Bogoliubov quasiparticles jk- 

l*(r)> = II ^i- k \o) , (8) 

a<k<-rr 

where |0) is the vacuum state of c-fermions. 
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B. Dynamics 

The dynamics of the Ising chain is conveniently de- 
scribed within the Heisenberg representation^. The 
Heisenberg representation of an operator 0(t) is defined 
by H {t) = U\t)G(t)U(i), and it evolves in time ac- 
cording to the equation of motion: 

ijO H (t) = U\t) ([0{t),H{t)] U(t) + i|o(t)) U(t) , 

(9) 

[•. •] denoting the commutator of two operators. 

Since the Ising Hamiltonian is quadratic in the c- 
fermions, the Heisenberg equations of motion for the op- 
erators Cf?(t) are linear and can be solved with a standard 
Bogoliubov-de Gennes approach. Their solution is writ- 
ten according to 



c k H (t)=u k (t)^-v*(t)^J 
c{ H (t)=v k (t) 1 ° k +ut(t) 7 °J 



(10) 



where 7° are the Bogoliubov operators that diagonalize 
7i(To) at the initial time, and the coefficients Uk(t),Vk(t) 
obey the time-dependent Bogoliubov-de Gennes equa- 
tions and are given by 



u k (t) 
v k {t) 



1 




Ufc(O) 


s «fc* J 


» fc (0) 




, »*(0) 


= <° 


-(C 


* ^ 






* ) ' 





(11) 



and 



(12) 



III. DENSITY OF KINKS 

The quantity we will analyze throughout this paper is 
the density of defects (or kinks), defined by: 



L ^ 



1 



(13) 



where the sum goes from j = 1 to N, for periodic bound- 
aries. Switching to the Heisenberg representation, one 
can compute the expectation value of the density of kinks 
p(t) at a certain time t. In the thermodynamic limit, this 
is defined by: 

p(t) = lim (*(r )|A^(i)|*(r )). (14) 

L — >oo 

In order to evaluate this expectation value, one has first 
to express A/" H (t) in terms of c-fermions, and then write 
them as combinations of quasiparticles 7^ that diagonal- 
ize the Hamiltonian before the quench, using Eqs. (|10p . 



In terms of the coefficients Uk(t) and Vk(t), the expecta- 
tion value is written as 

r dk 

Jo 2^ 

-i(u k (t)vl(t) - u* k (t)v k (t)) sink . (15) 



P(t) 



1- (\v k (t)\ 2 -\u k (t)\ 2 ) cos k 



Substituting the solution of the Bogoliubov equations for 
Uk{t) and Vk(t), Eq. (TTTj) , one obtains the following inte- 
gral expression for the density of kinks at zero tempera- 
ture: 



P(t) 



dk 
2^ 



1 - 2- 



1 + T n cos k 



(16) 



,r(r -r) sin 2 fc 



(4) 



n2e r 



1 



(2e£t) 



The generalization to finite temperatures T is quite 
straightforward; however, concerning the case T > 0, 
we are only interested in the equilibrium situation, to 
which we restrict ourselves. If the average is taken on 
the thermal state at temperature T, the density of kinks 
at equilibrium is defined by 

p£? = lim iTr[e-^ r ^l , (17) 

L^oo Z 

where e~^ H ^ /Z is the canonical ensemble of the system 
at a given temperature T = 1/(3, and Z = Tr[e _/3W ( r )] is 
the partition function. At finite temperatures, one has to 
take into account the possibility of having thermally ex- 
cited quasiparticles, according to the Fermi distribution 
function 

MT) = „„, a ] r , , ; , (18) 



exp(/3 e 



1 



where is the energy of the 7^ quasiparticle. The 



quasiparticle averages on the thermal state can thus 
be evaluated in the same way as on the ground state, 
using {il"tl) T = <7m7^)t = and {^ V ) T = (l - 
{lnlt)T)5n,v = n^T) S^.v, where (-)t denotes the ther- 
mal average and „ is the Kronecker delta. Proceeding 
in a way analogous to the zero-temperature case, one ob- 
tains this expression for p^?: 

dk 



eq 

Pt = 



2n 



Tl2^ 



l-(l-2n fc (T))(|^| 2 
i(l - 2n k (T)) («£(«£)* - (ul)*vl) sinfc 



cos A; 



(19) 



H and 



In the following, we will analyze in detail Eq. 
Eq. (TU). 



IV. BEHAVIOR AT EQUILIBRIUM 

Let us now start focusing on the expectation value of 
the density of kinks for the case in which system is time- 
independent, i.e. at equilibrium. In this situation the 
density of kinks is a time-independent quantity, and is 
given by the expression in Eq. (CG 
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A. Zero temperature 

At zero temperature the system is frozen in the ground 
state |W(r)), that is the vacuum of particles, Eq. ((8]), 
therefore the Fermi distribution function reduces to zero 
and the only non zero expectation value is (7^7^) = 1. 
In this case, the density of kinks can be straightforwardly 
obtained by putting Tq — T in Eq. (|16p , and is given by: 



eq 

Po — P T=a 



rn dk 
n 2vr 



1 - 



1 + T cos k 



Vi + r 2 + 2rcos/fc 



(20) 



This expression can be computed analytically, and is 
written in terms of complete elliptic integrals^, which 
arc defined as follows: 



ATA) ee I' du 

Jo V(l-« 2 )(1- 

n(c,A) 



Xu 2 ) 



du 



Jo (1 - cu 2 )^(l - u 2 )(l - Xu 2 ) 
Following Appendix [Al one gets: 

f 1-r 2 



(21) 
(22) 
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7T 

r 2 — 1 

Trr 



-n(r 2 ,r 2 ) 



for r < 1 
for r = 1 



The density of kinks at equilibrium and at zero temper- 
ature, as a function of the field strength T, is plotted 
in Fig. [1] with black solid curve. As one can see, po is a 
monotonic increasing function with T, and its first deriva- 
tive diverges at the critical point T c = 1. 



B. Finite temperature 

The density of kinks at equilibrium at finite tempera- 
tures is given by Eq. (|19[) . and can be reexpressed in the 
following way: 



eq 

Pt = 



dk 
2^ 



1 - 



1 + T cos k 



vT + T 2 + 2r cos k 



;i-2n*(T)) 



(24) 

Unlike the zero temperature case, it is not possible to 
derive a simple analytic expression for p^?, and one has 
to evaluate it numerically. Results are shown in Fig. [TJ 
where we plot p^? as a function of the transverse field 
strength T. As for the T — case, increases mono- 
tonically with T, towards the maximum value 1/2, that 
is reached in the limit r — > 00. The singularity in its first 
derivative at r c disappears for T > and is progressively 
smoothed out, as T increases. 




FIG. 1: (color online). Density of kinks p!j? for the thermal 
state at temperature T, as a function of the transverse field F. 
The black solid line denotes po , that is the equilibrium density 
of kinks in the ground state |\I/(r)). Inset: first derivative of 
plj? with respect to V. 



V. BEHAVIOR AFTER A QUENCH 

We now discuss in details the behavior of the density 
of kinks at T = after a quantum quench in the trans- 
verse field. As explained in Sec. HH the system is sup- 
posed to be initially in the ground state at r , and then 
evolved according to the quenched Hamiltonian, follow- 
ing Eq. @. Therefore, concerning the non-equilibrium 
situation, we specialize to the zero-temperature case, and 
denote with PQ(t) the density of kinks. This is explic- 
itly given by Eq. (flUl) . and can be written as the sum 
of a time-independent part and a time-dependent term, 
which decays in time and vanishes for asymptotically long 
times: 



(25) 



with p^\t) ™ 

separately discuss the behaviors of p [ g ' and Pq' (t) 



In the rest of this section we will 

(0) a „A J 1 )/ 



A. Asymptotic value 



Similar to the density of kinks at zero temperature, 
po, the asymptotic value p^Q can be expressed in terms 
of complete elliptic integrals A^(A) and fl(c, A), Eqs. (j2T|) 
and (12"2")) . We introduce the following notation for conve- 



<P! ee r 2 - 2r r + 1 , ^ 2 = r r 2 + r - 2r . (26) 

The expression of pg is then given separately for the 
following three cases with respect to the value of To, as 
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explained in Appendix [Bl 
(i) Case T > 1: 



,(°) 



2(rg-i)(r r-i) 
(r + r)(r§ -!)„/! 1 



1 i 

2 ~ 2^ 



K 



(A) 



(27) 



r 2 

1 



n t i 

' 1 

-2 _ ivr>2 _ i\2 



(r -r)(r 2 -i)(r 2 - 1) 2 /^!^ i 



r <pi^2 



n(( 



(ii) Case T < 1: 
(0) _ 1 _ 1 

P Q ~ 2 27T 



2r ( r -r ) (i-r 1 )^ ) 



(28) 



(r + r)(i-r 2 ) 



n(r^,r 2 ) 

2 \ / 1 t^2 \ 2 



(r -r)(i~r 2 )(i-r 2 ) 2 n ^ 2 ^ ^ 



Vl^2 

(iii) Case To = 1: 

i i + r i - r 2 



(0) 

pq 



2tt 



8vr\/r 



In 



i + r + 2Vr 

1 + T-2VT 



(29) 



In Fig. [5] we show the asymptotic value of the density 
of kinks after the quench, as a function of the transverse 
field r. We notice that, when the system after the quench 
is a paramagnet (r > 1), the density of kinks always 
decreases with decreasing T, for any value of the initial 
field r > 0. In that case, the shape of the curve is 
qualitatively the same as that for the ground state value 
Po, even if it has a smaller gradient. The situation is 
subtler for quenches ending in the ferromagnet, where the 
behavior of p~Q becomes non monotonic. As a matter of 

fact, decreasing T, Pq has a minimum at a certain value 
r < 1 and then it increases abruptly. 



B. Time dependent transient 

The relaxation of the density of kinks toward its 
asymptotic value can be unveiled by an asymptotic ex- 
pansion of the time-dependent part of Eq. (flTJ)) . pq (t), 
for large t. As explained in Appendix [Bl an explicit ex- 
pression for the asymptotic expansion can be given sepa- 
rately for the following cases, up to order 0(t~ 2 ) in time, 
(i) Case T ^ 1 and T ^ 1: 



(*) 



r-r 



(30) 



1 / , , 7T 

. cos 4 1 - r \t - - 



i 



VT+r(i + rv 
+o(t- 2 ). 



cos (4(i + r)t+ J) | 




FIG. 2: (color online). Asymptotic value p^Q of the density 
of kinks for t — » oo for a quench from To > 1 (left panel) and 
from To < 1 (right panel), as a function of the Y after the 
quench. The black dotted line shows the T — equilibrium 
density of kinks po . 



(ii) Case r ^ 1 and V = 1: 

i-r 



320r(i + r o ); 3 / 2 

(iii) Case r = 1 and T ^ 1: 



cos(WJ)+0(t- 2 ). (31) 



P { Q ] (t) = -i{i S in(4(l -F)l 



(32) 



+ 



v^(i-r) 



4 V / 2f(r 

o(r 2 ). 



T)t 3 / 2 



<s U(l + T)t- 



The decay of Pq ' (t) for large times is a power-law with 
oscillations. This feature does not depend on whether the 
quench is across the critical point or not. If the system 
before the quench is not critical, To ^ 1, the leading term 
decays as t~ 3 / 2 . The oscillatory part of the leading term 
consists of two frequencies, 4(1 + T) and 4|1 — T|; the 
second vanishes if the system is quenched to the critical 
point. On the other hand, if the system before the quench 
is critical, To = 1, the leading term decays as t~ l with an 
oscillating term of frequency 4|1 — T|. These frequencies 
of the oscillation come from the modes with zero group 
velocity, defined by de\/dk — 0. Indeed, such modes are 
found from Eq. ((BJ) to be fc = and k = 7T, the latter of 
which disappears when T = 1. As it is seen from Eq. (|16p . 
the frequency corresponding to the mode fc is given by 
2e£. This leads to the frequencies 4(1 + T) for k = and 
4|l-r| forfc = 7r. 

Figure [3] displays the behavior of pq ' (t) for various 
types of quenches, obtained by numerical integration of 
Eq. <|T~6|) (solid curves), as well as by the analytic asymp- 
totic expressions in Eqs. (|3"0)) - (f3"2"]) (dotted curves). As 
it can be seen from the main panel, if the system be- 
fore the quench is not critical, the power-law exponent of 
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FIG. 3: (color online). Absolute value of the time-dependent 
transient for the density of kinks after various types of 
quench. In the main panel we show two cases of quench from 
a non critical ground state, while in the inset we start from 
the critical point. Continuous curves are obtained from a 
numerical integration of Eq. (|16[) while dotted curves are plots 
of the analytic expressions (|30[l - (|32[ l. The blue dashed lines 



denote the power- law behaviors Pq (t) 



t 



-3/2 



(main panel) 



and Pn\t) ~ t 1 (inset), and are depicted as guidelines. 



FIG. 4: (color online). Asymptotic value of the density of 
kinks as a function of the effective temperature T e s. Dif- 
ferent colors stand for various values of F, as explained in 
the caption. Symbols denote the density of kinks Pq after 
a quench; at each value of To corresponds a different initial 
state, therefore a different effective temperature, as explained 
in the text. Straight lines denote the finite-temperature equi- 
librium values p^? =T . 



Hamiltonian: 



the decay is £~ 3 / 2 ; notice also that, in the case where the 
system is quenched outside criticality, the oscillations are 
superpositions of different frequencies; if it is quenched 
to the critical point, oscillations are more regular and 
consist of just one frequency. On the other hand, if the 
initial state is the ground state critical, the decay is t -1 , 
as shown in the inset. 



f 17 dk r i 
<W(r)) Trff = / — e T k n k {T eS ) + rc_ fc (T off ) - 1 (33) 
Jo 2tt L J 

where we used the representation of T~t(T) in terms of free 
fermionic Bogoliubov quasiparticles, Eq. ([7]), and n k {T) 
is given by Eq. (fT%|) . We therefore define the effective 
temperature by the implicit equation^: 



VI. NON-THERMAL BEHAVIOR OUTSIDE 
CRITICALITY 

We now concentrate on the asymptotic value pq of 
the density of kinks after a quench, and discuss the pos- 
sibility to track a behavior of this quantity in terms of an 
equilibrium situation at a finite temperature^. In par- 
ticular we would like to study if it is possible to define 
an effective temperature for the quenched system (out of 
equilibrium) in the most natural way so that, looking at 

p^q , the system behaves as if it was at equilibrium, at 
the same effective temperature. 

Due to the quench, the ground state |^(ro)) of a 
given Hamiltonian H.(Tq) becomes an excited state for 
the quenched Hamiltonian H.(T), therefore it has a posi- 
tive energy E = i(*(r )|7i(r)|*(r )), as compared to 
the ground energy (we take all the energies normalized 
per single site). An effective temperature for the system 
out of equilibrium can be defined by equating such energy 
Eq to that of a fictitious thermal state for the quenched 



E ee (H(T)) Tc!{ . (34) 

For a given pair of values (ro,r), this always admits a 
solution; we point out that, for a fixed T, there are two 
values of Tq corresponding to the same T e g , one for Tq < 
r and the other for To > T, except at T c g = 0, which 
coincides with the static case To = T. 

We are now ready to perform a quantitative compar- 
ison between the values of p^q after a quench and p e ^ {{ 
at equilibrium, where T g for the out-of-equilibrium sys- 
tem is obtained from Eq. ([34]) . This is done in Fig. [4] 
As one can see, outside criticality the two quantities are 
evidently not related and behave in different ways. On 
the other hand, at the critical point r c = 1 symbols (di- 
amonds) perfectly follow the solid line (green data). We 
interpret this as a thermal behavior, in the sense that the 
density of kinks after a quench to the critical point is uni- 
vocally determined by an effective temperature T e ff that 
depends only on the initial state energy after the quench, 
according to Eq. (|3~4"|) . Remarkably such behavior is not 
found for a non-critical dynamics, where fine details of 
the initial condition seem to be important. 
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The thermal behavior at criticality can be recovered 
analytically. Indeed, in that case the density of kinks at 
equilibrium with an effective temperature T e g is given, 
by substituting r = T c = 1 in Eq. (f2"4"|) : 



dk 
2^ 

(W(r)) Toff 



l-\e^[l-2n k (T cS )] 



En 

4 



(35) 



where e£ c 



2V2 + 2 cos k is the energy of the 7fc-fermion 
at criticality (r = T c = 1), the second equality follows 
from Eq. (|33p with rik{T) = n_fe(T), and the third from 
the definition of T e s in Eq. (j34"| . The energy E of the 
state | Vl/(ro)} is obtained directly from Eq. ([7|), by eval- 
uating (*(r )|7^7fc|*(r )), and is given by 



En = - 



* dk 2(1 + r )(l + cos/c) 



(36) 



On the other hand, the value of 



(0) 



for a quench at 



criticality reads, from Eqs. (|T6|) and lp5"]) . 



w dk 
2^ 



l + r cosfc (r -l)sin 2 fc 



rc^To 



which, after simple algebra, can be shown to reduce to 
1/2 + .Eo/4, hence obeying the rigorous equality 



- rt eq ' r<! 



(37) 



VII. CONCLUSIONS 

In this paper we have analyzed the behavior of the den- 
sity of kinks in the quantum Ising chain, both at equilib- 
rium and after a sudden quench in the transverse field. 
At equilibrium the density of kinks is monotonically in- 
creasing with the transverse field strength, and presents 
a divergence in its first derivative at the critical point. 
This singularity is smoothed out by finite temperature 
effects. In the quenched situation, it exhibits a tempo- 
ral decay down to a residual value which depends on the 
initial and the final values of the field. The finite-time 
transient is characterized by an oscillating power-law de- 
cay, whose rate depends on whether the initial state of 
the system is critical or not. 

We have also shown that, if the system is quenched 
towards the critical point, the density of kinks in the out- 
of-equilibrium situation can be described by equilibrium 
statistical mechanics, provided an effective temperature 
is defined according to the energy of the state after the 
quench. We say that this observable exhibits a thermal 
behavior only for a critical dynamics, while it does not 
thermalize for quenches to non-critical phases. These re- 
sults should be regarded in the light of the discussion in 
Ref^, where it was conjectured that only non-local op- 
erators in the fermion quasiparticles that diagonalize the 



model in the continuum limit exhibit a general thermal 
behavior. Indeed, we find non-thermal behavior for our 
density of kinks, which is a local operator in terms of 
fcrmions. However, the surprising thermal behavior at 
criticality should be better understood in a more general 
framework. 
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APPENDIX A: DENSITY OF KINKS IN THE 
GROUND STATE 

In this appendix we derive an analytic expression for 
the expectation value of the T = density of kinks po 
with respect to the system ground state |W(r)), in the 
thermodynamic limit. 

As explained in Sec. IIV A[ the density of kinks for the 
ground state is given by Eq. (|2"0"|) . Applying the variable 
transformation x = cos k we get 



Po 



dx 



1 + ra; 



1 1 



2 2vr J_ t Vl -x 2 v 7 ! + T 2 + 2Tx 



(Al) 



In order to evaluate the integral, we consider the follow- 
ing three cases with respect to T; 
i) Case T > 1. We perform the change of variable 

Tu+1 

x = TF ' 

u + 1 

so that the integral in r.h.s. of Eq. (|Alj) becomes 



dx 



1 + ra; 



i VI -x 2 VI + T 2 + 2Tx 

1 du r 2 - 1 

-1 r a/O 71 u 2 ){l-u 2 /T^ 



1 



1 -u 2 /T 2 



where complete elliptic integrals are defined by Eqs. ((2"T|) - 
Therefore we obtain: 



i r 2 -i 



P0 ~2 nF 



(A2) 



ii) Case T < 1. From the change of variable 

u + r 



Tu+ 1 



8 



it follows that: 

- 1 dx 1 + Tx 



i VI -x 2 VI + T 2 + 2Tx 



du- 



1 



1 -V 1 



i ^(l-^xi-r 2 ^ 2 ) i - r 2 w 2 
2(1 -r 2 )n(r 2 ,r 2 ). 



Therefore: 



1 1 -r 2 



Po = ■x 



n(r 2 ,r 2 ) 



2 7T 

iii) Case T — 1. The integral in Eq. (|A1|1 reduces to 
1 f 1 dx 



(A3) 



Therefore: 



V2J-1 VT 



1 1 

Z 7T 



(A4) 



APPENDIX B: ANALYTIC CALCULATION OF 
THE DENSITY OF KINKS AFTER A QUENCH 

We derive here analytic expressions for the density of 
kinks after a quench Pq{€). We consider both the time- 
independent part p*Q and the time-dependent transient 

pg } (*). 



1. Asymptotic value 

The time-independent part p'q of the density of kinks 
after a quench is written in integral form, starting from 
Eq. as 



JO) 



r dk( \ 
J 2^Y 



r cosfc r(r - r) sin 2 /c 



1 



(Bl) 



Applying a variable transformation by x — cos k and sub- 
stituting e\ = 2Vl + r 2 + 2r cos k, the integral / is ex- 
plicitly written as: 



1 



dx 



1 + T x 



2tt J_ t vr 



VI + r 2 + 2r a 



T{T -T){l-x 2 ) 



(1 + r 2 + 2r a; ) v / i + r 2 + 2r oa 



(B2) 



In order to evaluate the integral, we separately consider 
three cases with respect to IV 

i) Case To > 1. We apply a variable transformation 

r oM + i 

x = . 

u + T 



Then Eq. (|B2p is arranged into 

_ 1 r 1 du 1 



2^ J-i r 2 , V(i-« 2 )(i-« 2 /rg) 
r (r 2 -i)(r r-i) , (r + r)(r 2 -i) 



+ 



^ 2(i- u 2 /r 2 ) 
r (r -r)(r 2 -i)(r 2 -i) 2 - 



where we remark the notation defined by Eq. (|26p . Using 
the definition of the elliptic integrals, Eqs. (|21[) and 
we reduce it to 



1 r2(r 2 -i)(r r-i) r / 1 



(a) 



j = — — — 

, (ro + r)(r 2 -i)/ 1 1 



(B3) 



r 2 
1 



KrTf 2 ) 



+ (r -r)(r 2 -i)(r 2 -i) 2 ° n / yi 2 ^ j_ 

Hence, from Eq. (|B1[) . we get the analytic expression for 
pg } in Eq. M 

ii) Case < To < 1. We perform a variable transforma- 
tion by 

u + r 



Tqu + 1 

Using the new variable u, Eq. (|B2p is written as 



I = -J- 



2tt 7-i va-^Xi-rgu 2 ) 

r(r -r)(i-r 2 ) r + r i-r 2 



r 0V j 2 2r i-i> 2 

(r -r)(i-r 2 )(i-r 2 ) 2 ' 
2^1^2(1 - ifW/fi) 



The integrals above are expressed in terms of complete 
elliptic integrals to yield 



I = — 



+ 



1 ; 2r(r -r)(i-rg) 

2vr I r v?2 

n(r 2 ,r 2 ) 



(B4) 



(r + r)(i-r 2 ) 

To 

(r -r)(i-r 2 )(i-r 2 ) 2 TT /^ 2 ,2 ^ 



flf2 



n (©^)} 



Therefore, taking Eq. (|B 1 1) into account, one obtains the 

analytic expression for pq* in Eq. 1281 

iii) Case Tq = 1. Eq. (|B2[) is simplified into 



J=-i- 



rf.r 



1 r(i-r) i-x 

2tt VT^ I V^ + V2 1 + T 2 + 2Tx 



The first term in r.h.s. yields Regarding the 

second term, we apply a variable transformation: t = 



9 



r(i-r) 



x). Then integral is carried out and yields 

i + r + 2vT 



-4Vf + (1+T) In 



27rV2(2r) 3 / 2 
Summing the two terms, / is written as 



i + r-2Vr 



I = - 



i + T i - r 5 



2tt 



8vrVr 



hi 



i + r + 2Vr 
i + r-2Vr 



From Eq. (|B1[) . one obtains the analytic expression for 



Pq m Eq 



2. Time-dependent transient 



Following Eq. (|16p . the time-dependent part pg\t) of 
the density of kinks after a quench is written as 



(t)=8 



w dfc r(r - T) sin 2 k 



2tt 



cos(2e£i) 



Applying a variable transformation by u = 2 (r+i) 
equation is arranged into 



this 



P { Q>(t) 



(r -r)(i + r) 
47rVr^r 

1 dU y/(u 2 -Upjl-U 2 ) 



(B5) 



cos(t u), 



where we defined u = |1 - r|/(l + T), t' = 4(1 + r)t, 
and V = (r - r)(l - r r)/{r (l + r) 2 }. Since tp as a 
function of T has a derivative dip/dT = T(l — Tq)/Tq, 
it is monotonically increasing with To when < To < 1, 
while it decreases when To > 1. The maximum of ip 
for a given T is found at Tq = 1 and its value is (1 — 
r) 2 /(l + r) 2 = Uq. The variable ip can be either positive 
or negative, depending on T and To. 
i) Case To ^ 1 and T ^ 1. We introduce a number a that 
satisfies 



< a < minjito — He(y/tjj), 1 — uo}, 
where we remark that u Q — Re(y / V') is larger than zero, 



since < uq < 1 and < u 2 .. We define 



MO 



ha du v/(u 2 -u 2 )(l -u 2 ) 



1 



cos(t'u),(B6) 



du y (u 2 



<)(l-u 2 ) 



u +a 



V^ 2 



1> 



cos(t'u). (B7) 



We first consider the asymptotic behavior of I\{t') for 
large t' . Shifting the variable as u' = u — uq, the integral 
Ii(t') is arranged into 



h{t') = 



"2u (l - uD' 


1/2 




Jo 



du' COS [t'(u' + Uq)] 



(1 + Ji-)«'(1 - t^— )(1 + t^— ) 



(i + ^t) 2 (i + 



)(1 



1/2 

(B8) 



Since for < v! < a 



u' 


, and 


u' 









2uq ' 1 — Uq' 1 + Uq' Uq' 

are non-negative and less than unity, one can expand the 
integrand term inside the square brackets into a power 
series of the type [ • ]^ 2 = v / u 7 (l+aiu'+a 2 u' 2 + - ■ ■ )• Ap- 
plying this power series expansion and using asymptotic 
expansions of the Fresnel's type integrals 35 to Eq. (|B8|) . 
one obtains: 



h(t) 



"2u (l - u 2 )~ 


1/2 | 


_u 2 {u 2 -ip) _ 





sin[t'(«o + a )] 



a\a + a 2 a + • • • ) 



2t' 3 / 2 



t'u 



As a final result, after rcsumming the power series of a, 
the integral I\ (t) is evaluated as 



h(t) = 



(2u + a)a{l - (u + a) 2 } 



(u + a) 2 {{u + a) 2 - tp} 



2u (l - u§) 



1 /2 

sin[t'(u + a)] 



1/2 



2t' 3 / 2 



cos 



. u o( u o - i>) 

+0{t'- 2 ). (B9) 
We next consider h(t)- Transforming the variable by 



— 1 — u, Eq. (|B7|) is arranged into 
h{t) = 



[2(1 


-uir 


1/2 „\— Uq 


1 




Jo 



du' cos [t'(l- u') 
( 1 -lt^)( 1 -T^K( 1 -T) nl " 

Similarly to the evaluation of I\ (t) , one can define power 
series of u' for the term inside the square brackets: 
[■] x ^ 2 — v / u 7 (l + feiu' + & 2 'u' 2 + - ■ ■ ). Using asymptotic ex- 
pansions of the Fresnel's integrals^ 5 -, the integral is eval- 
uated asymptotically as 



h{t) 



\2{l-u 2 y 


1/2 | 


1-tp 





sin [i'(uo + a)] 



y/1 — Uo — a 



x(l + - u -a) + b 2 {l-u -a) 2 + ■■■) 



2t' 3 / 2 



t' + 



/-2' 



Finally, returning the power series into an original func- 
tion, one obtains 



h(t) = 



(2u + a)a{l - (u + a) 2 } 



-i 1 / 2 



(uo + a) 2 {(u a + a) 2 - tp} 
2(1 -u 2 ) 



l -v 
o(t'- 2 ). 



1/2 



2t> 3 / 2 



t' + 



in [t'(ua + a)] 



(BIO) 
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Looking at Eqs. (|B9[) and (|B10|) . one finds that the first 
terms in r.h.s. of both equations cancel in I\(t) + hit). 
After expressing uq, tp and t' in terms of T, Tq and t, 
from Eq. (|B5|) we arrive at the following expression: 



(*) 



r-r 



WV^rTt 3 / 2 



(Bll) 



1 



vH^rf|i-r | 
1 



■ cos 



4ji - r|t- - 



+o{r 2 ). 



■ cos 



4(i + r)t- 



ii) Case To ^ 1 and V = 1. Notice that uo = and 
ip=-(l- r ) 2 /4r < 0. We hereafter define ip 1 = -tp. 
P^q (t) is written as 



it) 



r -i 



h{t), 



where we have defined 



hit) 



1-u 2 , . 
em — , cos(8iu). 



o \Ju 2 + ip 1 



(B12) 



(B13) 



Changing the variable by it' = 1 — u, the integral is ar- 
ranged into 



hit) 



i + i>' Jo 



du' cos [8i(l 



-,1/2 



d n. vi nL ) 



For < u' < 1, one can expand the term under the square 
root in power series of u'\ [ ■ J 1 ^ 2 = \fu' '(1 + ciii' + C2U' 2 + 
• • • ) . This and the asymptotic expansions of Fresnel's 
integrals^ 5 , lead to: 



hit) = 



■ cos 



1" 

8i+ - 

4. 



+ 0(<~ 2 ) 



1 + V 32 v / 2t 3 / 2 
Writing ^' in terms of Tq, Eq. (|B12j) is then evaluated as 



(*) 



i-r„ 



32V^(i + r )t 3 / 2 



cos 



7f 

8t + — 

4 



+ o(t~ 2 ). 

(B14) 



iii) Case r = 1 and T ^ 1. Notice that = Uq. Eq. !fB"5|) 
is reduced to 



i-r 2 

47rVr 



hit), 



(B15) 



where we have defined 



hit) = / du 



1 , vr 



cos(i'u). (B16) 



U 



Changing the variable by u' = 1 — u, the integral is ar- 
ranged into 



l-ito 



hit) = V2 du' cos [t'il - u')] x 

Jo 



i - u' 



For < it' < 1 — uq, we perform a power series ex- 
pansion in v! as follows: [ • ] = 1 + diu' + d^u' 2 + ■ ■ • . 
Using this and the asymptotic expansion of Fresnel's in- 
tegrals 35 , h{t) is evaluated as follows. 



hit) = V2^- V1 U ° -sm(u t') 



x(l + di(l - u ) + <fa(l - uo) 2 + • • • ) 



2i' 3 / 2 



<+ 4 



0(t' 



The power series of (1 — uq) are summed to yield 

hit) = 



vv o; sin (w t ) - W-t-ttt^cos 



-o(t'- 2 ) 



2 t' 3 / 2 



*' + 



Writing uo in terms of T and substituting t' = 4(1 + T)t, 
Eq. (|B15[) is evaluated as 



PqW = -^{^in[4(l-r)t] 

0F(i-n „„„ 



(B17) 



4^(1 + r)i 3 / 2 
+o(t- 2 ). 



4(i + r)t +I 
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